Vertex-corrected tunneling inversion in superconductors: Pb 
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Abstract 

The McMillan-Rowell tunneling inversion program, which extracts the 
electron-phonon spectral function a 2 F(Q) and the Coulomb pseudopotential 
jx* from experimental tunneling data, is generalized to include the lowest-order 
vertex correction. We neglect the momentum dependence of the electron- 
phonon matrix elements, which is equivalent to using a local approximation. 
The perturbation theory is performed on the imaginary axis and then an ex- 
act analytic continuation is employed to produce the density of states on the 
real axis. Comparison is made with the experimental data for Pb. 
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The theory of low-temperature superconductors is one of the most accurate theories in 
condensed-matter physics. Agreement to better than one part in 10 4 is common between 
the tunneling density of states (DOS) measured experimentally and that calculated with 
an extracted electron-phonon spectral function a 2 F(Q) and Coulomb pseudopotential fi*. 
The reason why the agreement is so good is due to Migdal's theorem as formulated by 
Eliashberg || in the superconducting state, which says that there is a small parameter in 
the theory, namely the ratio of the phonon energy scale to the electronic energy scale, that 
guarantees the rapid convergence of the perturbative expansion. However in some recently 
discovered materials, that are believed to be electron-phonon superconductors, the ratio of 
the phonon to electronic energy scale is no longer as small. Two examples are Bai-^K^BiOs 
and the doped fullerenes 0] which have relatively high transition temperatures. This 
motivates the need for a theory that includes the effects of the so-called vertex corrections, 
which are neglected in the conventional Migdal-Eliashberg formalism. In this contribution, 
we present a generalization of the McMillan-Rowell |J tunneling-inversion program to in- 
clude the lowest-order effects of vertex corrections in the local approximation, and apply 
the new formalism to the low-temperature superconductor Pb to illustrate how vertex cor- 
rections can be incorporated into such an analysis, and to determine, quantitatively, the 
accuracy of the Migdal-Eliashberg formalism. 

Recent work on incorporating vertex corrections into the theory of superconductivity 
has illustrated some of their qualitative effects: If one adopts the conventional approx- 
imation of assuming a constant electronic DOS and neglecting the momentum dependence 
of the electron-phonon matrix elements, then the vertex corrections will suppress T c and the 
isotope coefficient a. Incorporation of either momentum dependence to the matrix elements, 
or nonconstant DOS can lead to enhancements to T c and a. Little is known about how large 
these effects can be in real materials, but they have been verified for model systems by 
comparing vertex-corrected theories || to the exact solution of electron-phonon models in 
the infinite-dimensional limit [|TD[] . 

Even a well-studied low-temperature superconductor such as Pb has a nagging inconsis- 
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tency between the experimentally extracted a 2 F and fi* and the bulk transition temperature 
T c . The extracted dimensionless electron-phonon coupling strength A satisfies A = 1.55, and 
the Coulomb pseudopotential is /i* = 0.131 ||. However, the Coulomb pseudopotential must 
be increased to fi* = 0.144 |Tl[| in order to produce the correct T c of 7.19K. It is possible 
that including the vertex corrections can explain this discrepancy. 

Furthermore, new computational strategies have been developed that greatly improve the 
accuracy of the Migdal-Eliashberg formalism, and allow for a straightforward generalization 
to the incorporation of the lowest-order vertex corrections: (i) use of the local approxima- 



tion for the many-body problem [fl^j ; (ii) performing exact analytic continuations from the 



imaginary axis to the real axis [13|; and (iii) incorporation of high-frequency resummation 



schemes |]T4| , p~5|| . We employ all three strategies in our computational methods. 

The electronic self-energy and the irreducible vertex function for superconducting order 
including the lowest-order vertex correction beyond the conventional model are both illus- 
trated in Figure 1. The wavy lines denote dressed phonon propagators, and the solid lines 
denote dressed electronic Green functions in the Nambu-Gorkov formalism. We make the 
conventional approximations |TT| : (i) neglect the angular dependence of the electron-phonon 
matrix elements and the phonon spectral functions, and evaluate them at the Fermi energy; 
(ii) neglect the energy dependence of the electronic density of states p(e) and evaluate it 
at the Fermi energy p(0); and (iii) treat the Coulomb interactions via a pseudopotential 
for the anomalous self-energy. When these approximations are invoked, one must solve a 
self-consistent perturbation theory in frequency-space only — the self-energy is replaced by 
a momentum-independent function that has been averaged over the Fermi surface. These 
self- consistent equations are identical in form to the equations one would derive in the local 
approximation, valid in the large-dimensional limit 

The perturbation theory is performed on the imaginary axis at the electronic Matsubara 
frequencies u n := irT(2n + 1). This allows for a proper treatment of the Coulomb pseudopo- 
tential, since the sharp cutoff lies on the imaginary, not the real axis . If we make the con- 



ventional definitions for the quasiparticle renormalization Z n := Z(iuj n ) = 1 — Im.[Yl(iu) n )]/u} n 



and for the superconducting gap A n := A(iu n ) = (p(iuj n )/Z(iuj n ), then the self-consistent 
equations are: 
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where A n := 2 / °° dVta 2 F (Vt)Vt / (Q 2 +An 2 T 2 n 2 ) is the dimensionless electron-phonon coupling, 
iV = \(*% — 1) is the cutoff for the summations [the frequency cutoff is chosen to be six 
times the maximal frequency in a 2 F(Q), or u c = 6uj ma x}, SZ n is defined below, and C is a 
Fermi-surface average for the vertex-correction terms. This average is defined by 
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with ep the Fermi energy and the band structure. In a free-electron model, with a k 
dispersion, the constant C assumes the form C = 4-, with n the number of free electrons 
per spin per unit cell J?]. Since the DOS at the Fermi level for a free-electron model is 
p(0) = the product Cp(0) is for a free-electron model. For Pb, we perform a 
scalar relativistic density-functional calculation of the band structure in the local-density 
approximation, and find p(0) = 2.5 x 10~ 4 states/sp in/me V, and C = 0.15, yielding Cp(0) = 
4 x 10 -5 states/spin/me V. This is the value we use for the numerical work. 

Finally, a high-frequency resummation scheme fl^UI^fl is employed, that calculates the 



perturbation-theory results relative to the exact results for the normal state, namely 

1 n 

Z n (normal) = 1 + — — [A + 2^A m ] , (4) 

^ n + 1 m=l 

by appending the normal state results to the perturbation theory illustrated in Figure 1. 
This is achieved by adding the contribution 
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5Z n : = Z n (normal) - 1 - ^ -^m s g n ( n ~ m ~ «) ' ( 5 ) 

m=—N ~ 

[where sgn(x) is 1 if x > and — 1 if x < 0] to the perturbation series. The term 5Z n 
is the normal-state contribution to the quasiparticle renormalization factor that is usually 
neglected when one introduces the cutoff N into the frequency summations. 

The difficult step in the tunneling-inversion program is to analytically continue the gap 
function to the real axis, in order to determine the tunneling DOS N(u) from 



N(u) = iV(0)Re 
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Recently, however, it was discovered that an exact analytic continuation could be performed 
by following the prescription of Baym and Mermin JTTj: formally perform the analytic 
continuation iu n —> u + iS, and then add a function that vanishes at each of the Matsubara 
frequencies, and is bounded in the upper half-plane except for simple poles located at just 
the positions necessary to cancel the poles introduced by the formal analytic continuation. 
Since the analytic continuation is unique, and since the final function can be shown to 
be analytic in the upper half-plane, this analytic-continuation procedure is exact. Such a 
scheme has already been implemented for the Migdal-Eliashberg theory [TB|], and it is a 



straightforward but tedious procedure to generalize these results to include the lowest-order 
vertex corrections. The end result is a formula for the quasiparticle renormalization, and 
the gap, on the real axis, that involves the data on the imaginary axis and integrals of the 
Green function evaluated in the upper half-plane. Because of the dependence on the Green 
function, these equations require a self-consistent solution on a real-axis grid (a step size of 
0.05 meV is chosen for the grid spacing). The final equations are cumbersome and will be 
shown elsewhere. 

Finally, we use this formalism to extract both a 2 F(Q) and fi* from the experimental 
data. We follow the original prescription of McMillan and Rowell ||: (i) Guess an initial 
value for a 2 F(Q). (ii) Adjust //* to reproduce the experimental superconducting gap at zero 
temperature Aq [which is defined from ReA(u;) = uo at uo = Aq]. An analytic continuation 



employing a Pade approximation is used to determine /i*, since A is determined to an 
accuracy of one part in 10 5 with such an approximation, (iii) Compute the functional 
derivative of the change in the tunneling DOS with respect to a change in a 2 F(Q). (iv) 
Determine the shift in a 2 F by solving the matrix equation 



This expression is discretized on the real axis, and a singular-value-decomposition is em- 
ployed to determine the shift Sa 2 F(Q), since there are small eigenvalues of the functional 
derivative matrix, which cause instabilities in updating a 2 F. (v) Determine the new a 2 F(Q) 
by adding a smoothed shift 5a 2 F(Q) to it, with a 2 F forced to behave quadratically in Q for 
Q < 0.5 meV. This procedure is iterated until convergence is reached. 

The results for the tunneling inversion for Pb for both the Migdal-Eliashberg theory, 
and the vertex-corrected theory are presented in Table 1. Various parameters are recorded 
including A, the "average" phonon frequency U[ n := exp[2 / °° dfl\ia(Q)a 2 F(fl)/flX}, the 
area A under a 2 F, fi*, A , T c , the maximum error, and the root-mean-square error of 
the fit. Note how a proper treatment of the cutoff for fi* and inclusion of the high-frequency 
resummation improves the calculated T c and that the vertex corrections modify A by the 
order of 1% even though the Migdal parameter [Cp(0)u max \] is on the order of 0.0007. A 
plot of the extracted a 2 F(Q) is given in Figure 2(a). The vertex-corrected fit is the solid line 
and the Migdal-Eliashberg fit is the dashed line. The difference between the two spectral 
functions a 2 F v (Q) — a 2 F ME (Q) is plotted in Figure 2(b). The vertex corrections produce 
slight enhancements to the transverse and longitudinal phonon peaks and they suppress the 
spectral weight in the region beyond the maximal bulk phonon frequency for Pb (about 9 
meV). Finally a comparison of the fitted tunneling DOS is compared to the experimental 
data || in Figure 3. The small dots are the experimental data points used in the fitting 
procedure; the larger dots (from data taken at a higher temperature) were not used in the 
fit. The solid line is the vertex-corrected DOS and the dotted line is the Migdal-Eliashberg 
DOS. Note that the two curves lie on top of each other for low energy, but deviate more at 
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higher energies. It is difficult to tell which fit is more accurate. 

We conclude with a discussion of the general properties of vertex corrections. In the 
conventional model, where the momentum dependence of the electron-phonon matrix ele- 
ments is neglected, and the electronic DOS is assumed to be constant, the vertex corrections 
will always cause a reduction to T c and the isotope coefficient a. A generalization of the 
formula for the Holstein model || shows that the ratio of the vertex-corrected T c to the 
Migdal-Eliashberg T c is 

f°° dQa 2 F{Q) + \ ndtta 2 F{tt) f°° dQ'a 2 F{Q')—^—] 
A \Jo X Jo Jo + U7 

(8) 

in the weak-coupling limit (if we set fi* = 0). This formula gives an order-of-magnitude 
estimate for when effects of vertex corrections should be important in a real material (for 
Pb the ratio is 0.9978). 

The other main effect of vertex corrections is that they will modify the higher-energy 
structure in a 2 F(Q), because they involve processes where two phonons scatter, so that 
the structures in a 2 F(Q) are changed at multiples of the lower-energy peaks. This effect is 
small in a low-temperature superconductor, but could be significant in the newly discovered 
superconductors. 

Finally, the vertex corrections will also modify the isotope coefficient. It is possible that 
effects of vertex corrections could be seen upon reexamination of isotope coefficient data 
on low-temperature superconductors. Work in this direction is currently in progress. We 
also plan on extracting a 2 F and /i* for other low-temperature superconductors to see if 
the effects of vertex corrections are important in other materials. We plan on applying the 
vertex-corrected tunneling inversion program to both Ba!_ x K x Bi03 and the doped fullerenes 
once accurate tunneling data for these materials becomes available. 
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TABLES 

TABLE I. Comparison of fitted tunneling inversion data for the Migdal-Eliashberg theory and 
the vertex-corrected theory. The vertex corrections modify A by the order of 1%. 



Theory 


A 


ijj hl (meV) 


A (meV) 


M* 


A (meV) T C {K) 


max. error 


r.m.s. error 


Migdal-Eliashberg 


1.542 


4.863 


4.029 


0.136 


1.400 7.23 


0.0004 


0.0001 


Vertex-Corrected 


1.561 


4.847 


4.070 


0.141 


1.400 7.22 


0.0007 


0.0001 
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FIGURES 



(a) 




FIG. 1. Feynman diagrams for (a) the electronic self-energy and (b) the irreducible vertex 
function in the superconducting channel. The first diagram in (a) and (b) is the Migdal-Eliashberg 
approximation, the remaining diagrams are the lowest-order vertex corrections. The solid lines 
denote dressed electron propagators and the wiggly lines denote dressed phonon propagators. 
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FIG. 2. (a) Electron-phonon spectral function, a 2 F(Q), extracted from the experimental 
tunneling data for Pb. The solid line is the vertex-corrected fit and the dashed line is the 
Migdal-Eliashberg fit. (b) Difference in extracted spectral functions a 2 Fy{Q) — oi 2 Fme{&)- Note 
the enhancement of the peaks and the suppression in the region where f2 is larger than the maxi- 
mum phonon frequency in bulk Pb. 
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FIG. 3. Comparison of the experimental tunneling DOS (solid dots) to the fitted DOS for Pb. 
The solid line is the vertex-corrected fit and the dotted line is the Migdal-Eliashberg fit. The small 
dots are the experimental data included in the fit, while the large dots are experimental data not 
included in the fit, but taken at a higher temperature (and hence less accurate). 
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